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Abstract. New high statistics, high resolution fixed target experiments producing 10^ 
- 10® fully reconstructed charm particles are allowing a detailed study of the charm sec- 
tor. Recent results on charm quark production from Fermilab fixed target experiments 
E-791, SELEX and FOCUS are presented. 



INTRODUCTION 

This review contains recent results from three Fermilab fixed target experiments 
dedicated to study charm physics: E791 and E871 (SELEX) of charm hadropro- 
duction and E831 (FOCUS) of charm photoproduction. 

The experiment E791 used a 500 GeV/c tt" beam incident on platinum and dia- 
mond target foils and took data in 1991-1992. A loose transverse energy trigger 
was used to record 2^^ interactions. Silicon microstrip detectors (6 in the beam and 
17 downstream of the target) provided precision track and vertex reconstruction. 
The precise location of production and decay vertices of long lived charm particles 
allowed to reconstruct over 2^ charm particles. 

SELEX used a 600 GeV hyperon beam of negative polarity to make a mixed beam 
of S~ and vr" in roughly equal numbers. The positive beam was composed of 92 
% of protons, and 8 % of 7r+. Interactions occurred in a segmented target of 5 
foils, 2 Cu and 3 C. The experiment was designed to study charm production in 
the forward hemisphere, with good mass and decay vertex resolution for charm 
momenta in the range 100-500 GeV/c. A major innovation was the use of online 
selection criteria to identify events that had evidence for a secondary vertex. Data 
taking finished in 1997. 

The FOCUS experiment used a photon beam with < >= 170 GeV on a Beryl- 
lium Oxide segmented target. Sucessor of Fermilab E687, FOCUS had upgrades 

1) This work is partially supported by CLAF/CNPq Brazil and CONACyT Mexico 

To appear, Proc. VII Mexican Workshop on Particles and Fields, MeridSiYuc. Mex., Nov. 1999 
janjos@lafex.cbpf.br 

'^'^ ecuautle@lafex.cbpf.br 



in vertexing, Cerenkov identification, electromagnetic calorimetry and muon iden- 
tification. The target segmentation and the use of additional silicon microstrip 
detectors after each pair of target segments were major improvements. The ex- 
periment took data in 1996-1997 and collected over 7 billion triggers of charm 
candidates. FOCUS has the largest charm sample available in the world, with over 
10^ fully reconstructed charmed particles. 

The data available from these three experiments will allow to make high precision 
charm physics, to study rare decays and to search for new physics. 



CHARM PRODUCTION AND CROSS SECTIONS FOR 

D MESONS 

Charm production is a combination of short and long range processes. Perturba- 
tive Quantum Chromodynamics (pQCD) can be used to calculate the parton cross 
section, the short range process. The long range process of parton hadronization 
has to be modeled from experimental data. The two processes occur at different 
time scales and should not affect each other, leading to factorization properties. 
The general expression for production of charm in hadroproductions is [1] 

a{PA,PB) =J2 dXAdXBf\^{xA, l^)ff{xB, l^)^i,j{as{l^),XAPA,XBPB) (1) 

where ff^{xA, A*) is the probability of finding parton type i inside the hadron A with 
momentum fraction xa, and /i is the scale at which the process occurs. The other 
hadron taking part in the interaction has f^{xB, A*)- The elementary parton cross 
section aij{as{lJ^),XAPA,XBPB) term is calculated by QCD, according to parton 
models. Although QCD is well defined theory, solutions to most problems are quite 
difficult. So, to calculate the parton cross section we usually do a perturbative series 
expansion in terms of the strong interaction coupling constant, ag, and calculate 
at leading order (LO), next-to- leading order (NLO) and so on. Nowadays most 
calculations include NLO. Some Feynman diagrams for hadropro duct ion of charm 
at LO and NLO are shown in Fig.l (more details can be found in reference [1]). 
Usually the differential cross sections is measured as a function of the scaled lon- 
gitudinal momentum (Feynman xp) da/dxp, and as a function of the transverse 
momentum squared da/dpf. A phenomenological parametrization of the double 
differential cross section often used to describe the experimental data is given by, 

where A , n and b are parameters used to fit the data. In fact kinematic consider- 
ations provide predictions for the power n at high xp and QCD calculations give 
average pt values comparable to the charm quark mass. 




FIGURE 1. Some of Feynman diagrams for charm hadroproduction: leading order (left) and 
Next-to-Leading Order (right side) 

E791 has measured the total forward cross section and the differential cross sections 
as functions oixp and from a sample of 88990 ±460 mesons. Fig. 2 shows the 
differential distributions and its comparison to theoretical predictions from QCD 
calculations at NLO [2] and to the Monte Carlo event generator, Pythia/Jetset [3] 
for cc and DD production. Hadron distributions are softer than cc distributions 
due to fragmentation. With suitable choice for the intrinsic transverse momentum 
kt of the incoming partons {kt — 1 GeV/c), the Peterson fragmentation function 
parameter (e = 0.01), and the charm quark mass = 1.5 GeV/c^, NLO D meson 
calculation provides a good match to the distribution and fair match to the 
xp distribution. The hadronization scheme implemented in Pythia/Jetset can be 
adjusted to fit the data. The large number of events make it possible to clearly 
observe a turnover point greater than zero {xc = 0.0131 ± 0.0038) in the xp distri- 
bution. The positive value provide evidence that the gluon distribution in the pion 
is harder than the gluon distribution in the nucleon [4]. The total forward neutral 
D meson cross section measured by E791 is a{xF > 0) = 15.4 ±2:3 /x6/nucleon. 
SELEX has also prehminary data on cross section. Fig. 2 (right) shows the xp 
distribution for + cc. from E791 [4] and the preliminary data from SELEX. Up 
to Xp = 0.5 data overlap and agree well. However SELEX has no strong evidence 
for rise at large as seen in the E791 data. The continuous line is a fit to the 
data points with the phenomenological parametrization in Eq. 2. 

Correlations in the production of Charm pairs 

Observation of both of the charm particles in a hadron-produced event can give 
additional information on the production process and allow to test QCD predic- 




FIGURE 2. production: from E791 compared with theoretical models as a function of xp 
(left), (central column) and compared with preliminary results from SELEX (right). 

tions since both longitudinal and transverse momenta of the charm particles and 
angular correlations are explicitly measured. 

In the simplest parton model the charm and anticharm particles are expected 
to be produced in opposite directions in the transverse plane. However if one 

assumes that the incoming partons have an intrinsic transverse momentum, kf, 
this will affect the transverse momentum of the heavy quark pair, its azymuthal 
correlations and the transverse distribution of a single quark. The partons entering 
the hard interaction are indeed supposed to have a non vanishing primordial kt, 
seen as a nonpcrturbative Fermi motion of partons inside the incoming hadrons. 
Typical values of kt should thus be 300-400 McV. However it has been noted [5] that 
much higher values of kt are required, at or above 1 GeV, to reproduce charm data. 
This could be an indications of the importance of next to leading and higher order 
effects, by which emitted gluons would further modify the nearly back-to-back [6] 
production of the final charm hadrons. 

E791 has measured correlations between D and D mesons from 791 ±44 fully 
reconstructed charm meson pairs [7]. The main variables used to describe charm 
pair correlations are the beam direction, pt, and either the rapidity y or the Feyn- 
man scaling variable xp, and also the azymuthal distribution (f). The same way 
the difference and sum between these variables for two charm mesons {D, D) Ax^, 
TiXp and so on. 

The measured distributions are compared to predictions of the fully differen- 
tial NLO calculation for cc production [2], as well as to predictions from the 

Pythia/Jetset Monte Carlo event generator for cc and DD production. Default 
parameter have been used for the theoretical models. 



For the single charm distributions shown in Fig. 3, we observe that for the longi- 
tudinal momentum distributions xp and y the experimental results and theoretical 
predictions do not agree. In this comparison the experimental distributions are 
most similar to the NLO and Pythia/Jetset cc distributions, but are narrower than 
all three theoretical predictions. The experimental pi distribution agrees quite 
well with all three theoretical distributions. As expected, both the theoretical and 
experimental distributions are consistent with being flat. 

The experimental and theoretical longitudinal distributions for pairs Ax^,' , Sx^, 
Ay and are shown in Fig. 4. As with the single-charm distributions, the exper- 
imental results are much closer to the two cc predictions than to the Pythia/Jetset 
DD predictions, but narrower than all three theoretical predictions. For the trans- 
verse distributions for charm pairs any observed discrepancy between data and 
theory must derive from the theory modeling the correlations between the trans- 
verse momentum of the two D mesons pt,D and ^ because the single charm pi 
and experimental distributions agree well with the theory. The A0 distribution 
shows clear evidence of correlations (more details can be found in reference [7]). 




FIGURE 3. Single-charm distributions for XF^ y, pi and cf): weighted data (•) NLO QCD 

prediction ( ); Pythia/Jetset charm quark prediction ( ); and Pythia/Jetset D 

meson prediction ( ). All distributions are obtained by summing the charm and anticharm 

distributions from charm-pair events. 




FIGURE 4. Charm-pair Axf, T,xf, Ay, Sy, A^, S0, Ap^ and distributions; weighted 

data (•); NLO QCD prediction ( ); Pythia/Jetset charm quark prediction ( ); and 

Pythia/Jetset D meson prediction ( ). 



HADRONIZATION AND PARTICLE - ANTIPARTICLE 

ASYMMETRIES 

The production of a charm hadron can be subdivided in two steps: the production 
of a cc pair followed by the hadronization of these quarks. In perturbative QCD, 

that describes the cc production, the xp spectra of produced charm/anticharm 
quarks are identical to leading order and the effects of higher orders are very small 



in this respect. Therefore any asymmetry between charm and anticharm hadrons is 
a simple measure of nonperturbative effects coming from the hadronization process. 
Particle - antiparticle asymmetries can be quantified by means of the asymmetry 
parameter 

. N-N 

A = ^ 3 

N + N ^ ' 

where N (N) is the number of produced particles (antiparticles). This parameter 
is usually measured as a function oi xp and p^. 

Several experiments have reported an enhancement in the production rate of 
charm particles having valence quarks in common with the incident particles, rela- 
tive to charge conjugate particles which have fewer or no common valence quarks. 
This effect is known as leading particle effect. Measurements of the asymmetry 
parameter A can put in evidence leading particle effects, as well as other effects 
like associated production of a meson and a baryon. 

From the theoretical point of view, models which can account for the presence of 
leading particle effects use some kind of non-perturbative mechanism for hadroniza- 
tion, in addition to the perturbative production of charm quarks. Two examples 
are: 

String fragmentation [8]: in this case the parton of the hard interaction and the 
beam remnants are connected by a string which reflect the confining color field. 
Successive breaking of the color flux tube stretched between a cluster, when it is 
kinematic possible, will create light quark- antiquark and hadrons will be produced. 

Intrinsic charm model [9]: a virtual cc pair pops from the sea of the beam par- 
ticle. The cc pair coalesce with the neighbor valence quarks due to their similar 
rapidity. This mechanism favors the production of charm particles with valence 
quarks in common with the beam particle at high xp and low pt region. A similar 
argument can be drawn with respect to the target. 



New results on the asymmetry parameter A and evidence for leading particle 
effects in both meson and baryon production are available from experiments E791, 
SELEX and FOCUS. 

E791 has measured recently the asymmetries of Df mesons [12]. Fig. 5 shows the 
Df asymmetry as a function oixp and p^, compared with previous results from 
the same experiment [13]. 

Prehminary results of D^(D'^) as well as asymmetries as a function oi xp 
presented by SELEX [10] are shown in Fig. 6 (left), for different incident beam 
particles (vr", S^, p). The 7r~ data is compared to asymmetries from E791. 
The asymmetry at xp > 0.4 does not rise steeply with xp as previously reported 
by E791, Fig. 6 (right). 

For baryons there are also preliminary results for the asymmetry parameter. Fig. 7 

(left) shows the E791 results for the asymmetries as function of xp and p^, 
compared with predictions from Pythia/Jetset (full lines) [11]. The results show 
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FIGURE 5. Asymmetries of Ds as a function oi xp and , compared with . 




FIGURE 6. Asymmetries of D mesons from SELEX (preliminary results) and E791. The first 
and second columns are D^{D^) and asymmetries from SELEX. The right figure shows the 
comparison of asymmetries from SELEX and E791. 

a uniform positive asymmetry of 12.7 ± 3.4% over the studied kinematical range 
but do not exclude a rise in the xp < ^ region as predicted by Pythia/Jetset. For 
xp > Q the observed asymmetry does not agree with Pythia/Jetset predictions. 
SELEX has also measured the A J asymmetry as a function of Xi;' for different 



incident beam particles (vr^, S^, p) Fig. 7 (right) [10]. The asymmetry is clearly 
larger for the baryon beams than for the vr beam. For the protons the only region 
in which there is Ac production is at very small xp- Their preliminary results for 
the 7i~ beam are compatible with those from E791, Fig. 7 (left). 




FIGURE 7. Ac asymmetries as function oi and xp from E791 (left). Ac asymmetries (pre- 
liminary results) as function of xp for different beams {n~ , E^, p) from SELEX (right). 

FOCUS has also preliminary results on baryon asymmetries. In Fig.8 we see the 
Ac asymmetry as functions of pi, pf and xp obtained from a sample of about 16,000 
AcS, compared with Pythia/Jetset Predictions. As FOCUS has a photon beam, no 
leading particle effect is expected in the xp > region. In this case the positive 
asymmetry observed in all the xp range can be an indication of charm baryon and 
charm meson associated production, favouring a positive asymmetry. 
The high statistic Ac sample from FOCUS allowed to obtain about 600 Sc Acvr. 
It is interesting to compare the asymmetry for charm particles with different light 
quark content. We present in Fig. 9 (right) very preliminary results from FOCUS 
comparing the E++(m'uc) and T,^{ddc) to the A'^{udc) total asymmetry. 
In Fig. 9 (left) we show the comparison between the A° and the Ac asymmetries as 
a function of from E791. Their similarity suggests that the ud diquark shared 
by the produced A° (A+) and nucleons in the target should play an important role 
in the measured asymmetry in the xp < region. However, one expects that Ac 
asymmetry grows more slowly than the A° asymmetry due to its higher mass. 

Leading particle effects were also seen by E791 in hyperon production. Prelimi- 
nary results on A, S and Q asymmetries as a function of xp and p^ are shown in 
Fig. 10 in comparison with predictions from Pythia/Jetset. The range of xp covered 
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FIGURE 8. A+ asymmetry from FOCUS as a function of pi, p| and xp (preliminary results), 
compared with Pythia/Jetset predictions 
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FIGURE 9. Comparison between A° and Ac asymmetries from E791 (left). Ac, S++ and S° 
total asimmetries (preliminary) from FOCUS (right). 



allowed the first simultaneous study of the asymmetry in both the negative and 
positive xp regions. We can clearly see leading particle effects associated with the 
target or with the beam particles which qualitatively agree with expectations from 
recombination models (sec table 1) [14]. It is interesting to observe, as expected, 
the crossover of the S asymmetry with respect to the A asymmetry a.t xp — 0. 
The positive asymmetry measured in regions xp > for the A{udc) and for the 
Q{sss) suggest the associated production of a hyperon and a kaon due to the higher 
energy threshold imposed by baryon number conservation for the production of an 
anti-hyperon. Pythia/Jetset does not reproduce the data. 



TABLE 1. Hyperon asymmetries predictions 



X/ < 2;/ > 

target (uud or ddu) beam tt" {ud) 



A^{uds) double leading leading 

A^{uds) non leading leading 

E~{dss) leading leading 

S^(Jss) non leading non leading 

n~{sss) non leading non leading 

fl^ (sss) non leading non leading 

Recomb. Models'^: Aa > A-e > Ag As > Aa Aq 



^ Presented in Second Latin American Symposium [14] 
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FIGURE 10. Hyperon asymmetries as a function of (top) and (bottom) from E791 
(preliminary results). The asymmetry for Xp (p^) is integrated over all the p"^ (xp) range of the 
data set. The right column show the predictions of Pythia/Jetset both as a function of xf and 
p"^ (preliminary results). 



DOUBLE CABIBBO SUPPRESSED DECAYS 

The Cabibbo suppressed charm decays can provide useful insights into the weak 
interaction mechanism for nonleptonic decays. The D'^ — > K'^n'n'^ signal ob- 
tained from 100 % of FOCUS data set consist of ~ 300 events and is at least a 

factor of five larger than two previous observations by E687 and E791 ( E791 ob- 
served 59 ± 13 events [15]). The preliminary branching ratio relative to K~T:~^Tr~ 



is (0.72 ± 0.09) %, completely consistent with the world average of (0.68 ± 0.15)% 
and the E791 values of (0.77 ± 0.17 ± 0.08)%. We note that this is ~ 3tan^{0c), 
which is roughly the ratio of the / lifetime, indicating that the destructive 
Pauli interference present in the Cabibbo Favored decay is absent in the doubly 
Cabibbo suppressed (DCS) mode. 

— > K~K'^K'^ is an interesting DCS decay, which cannot even occur through a 
spectator diagram. FOCUS has the first observation of this mode, and reports a 
preliminary result for the Branching Ratio relative to K^n^n^ of (0.14 ±0.02) %. 
Several groups have reported observations of a — > 0/^+ signal, however FOCUS 
did not find evidence for such decay [16]. 

SELEX announced the first observation of a Cabibbo suppressed decay of a charm 
baryon through the decay pK~Ti^ [17]. Fig. 11 shows the signal of 157 ± 22 

events reported by SELEX and simple spectator diagrams with external W emis- 
sion for S+ decaying into a Cabibbo allowed and into a single Cabibbo suppressed 
(SCS) mode. The other Cabibbo allowed S~ mode interchanges s and d quarks 
lines and produces a dd pair from the vacuum instead of a du pair. FOCUS has 
also observed the same SCS decay, reporting a signal of 86 ± 21 events from about 
70 % of their data. 




FIGURE 11. Simple spectator diagrams (a) and (b) and preliminary signal of pK tt"*" 

from SELEX. 



E791 has pubhshed [18] results on the singly Cabibbo suppressed decay, D'^ — > 
K~ K~^TT~TT'^ . A coherent amplitude analysis of the resonant substructure was 
used to extract decay fractions. Significant phase angles among different modes 



TABLE 2. Branching fractions relative to 
D° K--K+n-'K+ from E791. 



Mode 


Branching Fraction 




(2.0 ±0.9 ±0.8) % 




(0.9 ±0.4 ±0.5) % 




< 2.0 % (90 % CL) 


K*OK+TT- + K*OK-Tr+ 


< 2.0 % (90 % CL) 



indicate very strong interference. The measured branching fractions relative to 
'tt'^tt tt"*" are presented in the table 2. 



Hadronic charm decays, Dalitz plot Analysis 

With the advent of high statistics experiments, charm meson decay have become 
a new way to study light meson spectroscopy. The amplitude analysis performed 
on Dalitz plots gives insight into the decay dynamics, providing direct information 
about intermediate resonances and relative decay fractions, and allowing to study 
final state interactions coming from the interference of the amplitudes describing 
competing resonant channels. 

E791 has preliminary results on the decay of D'^ and Df mesons in three pions. 
A clear signal with 1240 ± 51 and 858 ± 49 was obtained after apply- 
ing selection criteria aimed at identifying a clearly separated Svr vertex and after 
carefully estimating the backgrounds coming from possible reflections and three 
pion combinations. The branching ratios were normalized to D'^ — > K'tt'^tt'^ 
(34, 790 ± 232 events) and to 07r+ (1038 ± 44 events) respectively Efficien- 

cies were obtained from a full Monte Carlo simulation. The branching ratio of the 
D+ 7r+7r-7r+ relative to D+ K-ti+ti+ obtained was 0.0329 ± 0.0015l[i:[||j^i. 
Similarly the branching ratio of the — > 7r"'"7r~7r~ relative to D'^ — > 07r+ was 
0.247 ±0.028ig:g}^. 



Dalitz plot results from E791 

Among the advantages of using charm meson decays to study light 1=J=0 
states is the fact that, unlike hadron-hadron scattering, in the decays of D 
mesons the initial state is always — 0~, limiting the number of possible fi- 
nal states. The decay Tr^Tr+Tr"*" is Cabibbo-favored without a strange me- 
son in the final state. It can proceed via spectator amplitudes producing inter- 
mediate resonant states with hidden strangeness like the /o(980) or it can pro- 
ceed via W-annihilation amplitudes producing intermediate resonant states with 
no strangeness. The decays hke Df — > p^TT"*" and the non-resonant Df — > 7r~7r+7r+ 
would proceed via VT-annihilation mechanism. It would also be responsible for the 
decay —>■ /o(1370)7r"'", if the /o(1370) resonance is at least partially a uu + dd 




FIGURE 12. The tt+tt 7r+ Dalitz plot and its projections ( preliminary results) from 

E791. 



TABLE 3. Preliminary results of fo mass and width of Df /o(980)7r+ and 
D+ /o(1370)7r+ from E791. 



/o(980) 




/o(1370) 




Mass (MeV/c2) 


width (MeV/c2) 


Mass (MeV/c2) 


width(MeV/c2) 


E791 978 ± 4 


44 ±5 


1440 ± 19 


165 ±29 


PDG 978 ± 10 


40 - 100 


1200 - 1500 


200 - 500 



state as predicted by the simple quark model. The scale of the ly-annihilation com- 
pared to the ly-radiation amplitude would be indicated by the relative contribution 
of these channels to the tt^ti^ti^ final state [19]. 

The Dalitz plot of -D^ ti^h^h^ and the vrvr mass projections [20] are shown 
in Fig. 12. 

The /o(980)7r^ mode is the dominant one, accounting for nearly half of the 
D'^ 7r^7r+7r+ decay width. The /o(980)7r='= is often supposed to have a large 
ss component, indicating a large spectator amplitude in this decay. Significant 
contributions of /o(1370)7r+ and /2(1270)7r+ components were also found. The 
contribution of p^{770)n~^ and p''(1450)7r+ components correspond to about 10 % 
of the 7r^7r+7r+ width. This could indicate either contribution from the annihila- 
tion diagram or from inelastic final state interactions. No significant non-resonant 
component was found. 

Preliminary /q masses and widths [20] from E791 and PDG are presented in the 
table 3. 

FOCUS also has preliminary results of this decay mode. Their preliminary Dalitz 
plot shows the 7r^7r+7r+ based on a very clean signal of ~ 1300 events 

reconstructed from 100 % of their data. The preliminary results indicate a negligible 
contribution from the p suggesting negligible Weak Annihilation contribution [21]. 



FIGURE 13. £'+7r+7r~7r+ Dalitz plot and its projections. The central plot is without a aw 
state and the last plot include a cttt state in the fit (preliminary results). The shaded area is the 
background distribution. 

Dalitz plot results and evidence for a light scalar resonance 

The Dalitz plot of the single Cabibbo-suppressed decay 7r~7r+7r+ from 

E791 data is shown in Fig. 13 (left). A coherent amplitude analysis was used to 
determine the structure of its density distribution. The fit including a non resonant 
amplitude and amplitudes for decaying to a 7r+ and any of the five established 
TT+TT- resonances p°(770), /o(980), /2(1270), /o(1370), and p°(1450) is shown in Fig. 
13 (central). This fit is poor in the low vr+vr" region and has several unsatisfactory 
features [22]: the NR channel dominates, different from the decay, and the 
p°(1450)7r"'" is more significant than the p°(770)7r+ state. 

It was found that allowing an additional scalar state, with mass and width un- 
constrained improves the fit substantially. The mass of the resonance found by 
this fit is 486;26 MeV/c^ and the width 35i;43 MeV/c^. Referring to this Tr+Tr" 
resonance as the (t(500), it was found that —>■ cr(500)7r"'" accounts for about 
half of the total decay rate, non- resonant decay was very small and the p°(1450)7r+ 
fraction was much less than p°(770)7r"^. Preliminary results of the fit with this state 
are shown in Fig. 13 (right side). 

Theoretically, light scalar and isoscalar resonances are predicted in models for 
spontaneous breaking of chiral symmetry, like the a linear model [23]. These scalar 
particles have important consequences for the quark model, for understanding low 
energy vrvr interactions and also for understanding the AI = 1/2 rule. 



Multidimensional analysis of A+ pK 7r+ from E791 

E791 has reported recently the first amplitude analysis of the decay of a charm 
baryon [24,16]. The study of charm baryon decays can give information regarding 




FIGURE 14. A+ multidimensional analysis from E791: left plots shows projections of two body 
masses and angles. Solid line histograms are data in the range 2265 < MpK-^+ < 2315 MeV/c^. 
Dashed lines represent the fit in the same range. The right figure shows polarization Py+ vs. Pf. 



the relative importance of spectator and exchange amphtudes. Exchange amph- 
tudes are small in charm meson decays because of helicity suppression. However 
in charm baryon decays this effect should not inhibit exchange amplitudes due 
to be three body nature of the interaction. The spectator and W-exchange dia- 
grams can contribute to pK*'^(890), A(1520)7r"'" or pK~7r~^ modes. However for the 
A++(1232)fC^ the W^-exchange is the only diagram possible. 
The charm baryon can be produced polarized and its decay products carry spin. 
These extra quantum numbers require five kinematic variables for a complete de- 
scription of the decay, and instead of the conventional two dimensional Dalitz plot 
analysis, a five-dimensional amplitude analysis is required. 

A sample of 946 ± 38 A+ pK^n^ reconstructed decays was used by E791 to 
determine relative strengths and phases of resonances in the final state as well as 
the Ac production polarization. The fit projections and the polarization as function 
of pt are shown in Fig. 14. The resonant fractions for A+ pK~7r'^ are shown in 
table 4. 

The A"'"+(1232)ii'~ and the A(1520)7r"'" decay modes are seen as statistically signif- 
icant for the first time. The observation of a substantial A++(1232)ir~ component 
provides strong evidence for the W-exchange amplitude in charm baryon decays. 
It was also observed an increasingly negative polarization for the Ac as a function 
of Pt. 



TABLE 4. Resonant fractions for 
A+ decay modes. 



Decay Mode % of pK 7r+ 



pK*"{890) 
A++{1232)K- 
A(1520)7r+ 
non-resonant 



19.5 ±2.6 ±1.8 
18.0 ±2.9 ±2.9 
7.7 ± 1.8 ± 1.1 
54.8 ±5.5 ±3.5 



RARE AND FORBIDDEN DECAYS 

In the charm sector the rare and forbidden dilepton decay modes can be classified 
mainly into three categories (example of Feynman diagrams are shown in Fig. 15): 

1) Flavor Changing Neutral Current decays (FCNC) such as ^ and 

h+l+l-. 

2) Lepton Family Number Violating decays (LFNV) such as h'^lil2 and 
D° — > lfl2 where the leptons are from different generations. 

3) Lepton number Violating decays (LNV) such as D'^ — > h'l'^l'^ where the leptons 
are of the same generation but have the same sign. 

Where h stands for vr, K and / for e, /i. 

The first decay modes (FCNC) are rare, that means a process suppressed via the 
GIM mechanism which proceeds via an internal quarks loop in the Standard Model 
[25]. The FCNC decay mode l^l" can proceed via a W box diagram and 

the theoretical estimates [26] for the branching fraction are of the order of ~ 10~^^. 
The predictions for the other FCNC decay modes, h^l^lr . are considerably 

larger, of the order of ~ 10~^ . These decay modes can proceed via penguin dia- 
grams [25] and from long distance effects. 

The decays modes LFNV and LNV are strictly forbidden in the Standard Model as 

they do not conserve lepton number. However, some theoretical extensions of the 
Standard Model predict lepton number violation [27], and then the observation of 
a signal in these modes would be evidence for new physics beyond the SM. 

E791 has recently published [28] (see figure 15) a set of new limits in rare or 
forbidden decay modes that improve the PDG98 numbers by a factor of 10. They 
searched for 24 different rare and forbidden decay modes and have found no evi- 
dence for them. They therefore presented upper limits on their branching fractions. 
Fourteen of their limits represent a significant improvement over previous results 
and eight are presented for the first time. 

For this study E791 used a blind analysis technique. The mass region where 
the signal is expected is masked throughout the analysis. Selection criteria are op- 
timised by studying signal events generated by Monte Carlo simulation and back- 
ground events obtained from data in mass windows above and below the signal 
region. The criteria were chosen to maximize the ratio Ng/ ^J\Nb), where A^^ and 
Nb are the numbers of signal and background events, respectively. Only after this 




FIGURE 15. Rare decays limits from E791 (top) and some Feynman diagrams for FCNC, 
LFNV and LNV decay modes (bottom). 



procedure were events within the signal window unmasked. This bhnd technique 
is used so that the presence or absence of signal does not bias the choice of the 
selection criteria. 



FOCUS is also looking for rare and forbidden decays using the same technique. 
Some examples of rare decays are presented in tabic 5 where we compare the ex- 
pected sensitivity from FOCUS to E791 results and PDG values. 



TABLE 5. Rare decays limits from FOCUS (preliminary), E791 and PDG98. 



Decay mode 


FOCUS expected sensitivity 


£791^^ 


PDG 98 




(45 % of the data) 


90 % C.L. limit 


90 % C.L limit 


D+ K+n+ii- 


8.1 xlO-^ 


4.4 xlO-5 


9.7 xlG-5 


D+ K-ij,+ij,+ 


12.1 xlO-6 




1.2 xlO"'' 


D+ -K+jj+n- 


7.8 xlO-6 


1.5 xlO"^ 


1.8 xlO"^ 


D+ 'K-fJ,+ fJ,+ 


7.1 xlO-6 


1.7 xlO"^ 


8.7 xlO"^ 


D+ fJ.~fJ-+^+ 


4.4 xlO"'^ 







^ Published [28] 



SEMILEPTONIC DECAYS, FORM-FACTORS 

The weak decays of hadrons containing heavy quarks are influenced by strong in- 
teraction effects. Semileptonic charm decays such as — > K C^Ue, Df — > (pl^ui 
are an especially clean way to study these effects because the leptonic and hadronic 
currents completely factorize in the decay amplitude ,A, as we can see in the Eq. 
4, where Gp is the Fermi coupling constant for the weak interaction and Vcs is 
the CKM matrix element. L'^ (Eq. 5) and i?^ (Eq. 6) represent the leptonic and 
hadronic currents 



A{D+ ^ K*'e+iy,) = ^V^sL^H, (4) 



L^ = «e7''(l-75K (5) 



= {mo + mK*)Ai{q^)emu \ ■ Pd){pd +Pk)h - 

rriD + tjik* 



e ■ Pd) {pd - Pk)^ - i ■ ^ixupae^'PoPK (6) 

TJiD + rriK* mD + rriK* 

With a vector meson in the final state, there are four form factors, l^(q'^), Ai(g^), 
^2(5^) and A^lq'^), which are functions of the Lorentz-invariant momentum trans- 
fer squared q^, the square of the invariant mass of the virtual W [29]. The dif- 
ferential decay rate for D'^ — > with — > K'tt'^ is a quadratic ho- 
mogeneous function of the four form factors. Unfortunately, the limited size of 
current data samples does not allow precise measurement of the g^-dependence of 
the form factors; we thus assume the dependence to be given by the nearest-pole 
dominance model: F{q'^) — F(0)/(1 — q^/rn?poie) where rUpoie — rny — 2.1 GeV/c^ 
for the vector form factor V (which correspond to — l"*" state, -D*i), and 
TTipoie = "TTT-A = 2.5 GeV/c^ for the three axial- vector form factors A [30] (corre- 
sponding to 1~ state, D*). 



The third form factor A^lq"^), which is unobservable in the hmit of vanishing lep- 
ton mass, probes the spin-0 component of the off-shell W. Additional spin-flip 
amplitudes, suppressed by an overall factor of m| / when compared with spin no- 
flip amplitudes, contribute to the differential decay rate. Because v4i(g^) appears 
among the coefficients of every term in the differential decay rate, we can factor 
out Ai{0) and measure the ratios: 

ry = V{0)/Ai{0), r2 = A2{0)/Ai{0) and rg = As{0)/Ai{0). The values of these 
ratios can be extracted without any assumption about the total decay rate or the 
weak mixing matrix element Vcs- 
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FIGURE 16. Top-left figure present the ratios r2 and ry for Z?+ K*^lvi, and top-right 
compare them to theoretical models. Bottom- left show ratios 7'2 and ry for _D+ — s- (filvi, and 
bottom-right is a comparison with theoretical models. 



We report E791 measurements of the form factor ratios and r2 for the muon 
channel [31] and combined with electron channel [32] (see Fig. 16). This is the first 
set of measurements in both muon and electron channels from a single experiment. 
We also report the first measurement of = A3{0)/Ai{0), which is unobservable 
in the limit of vanishing charged lepton mass. 



The measurements of the form factor ratios for K jji'^y^ presented here 

and for the similar decay channel D"*" — > K^e^v^ [33] follow the same analysis 
procedure except for the charged lepton identification. Both results in the electron 
and muon channels are consistent within errors, supporting the assumption that 
strong interaction effects incorporated in the values of form factor ratios do not 
depend on the particular leptonic decay. 

The combined results of electronic and muonic decay modes produce ry = 1.87 ± 
0.08 ± 0.07 and r2 = 0.73 ± 0.06 ± 0.08. The combination of all systematic errors 
is ultimately close to that which one would obtain assuming all the errors are 
uncorrelated. The third form factor ratio = 0.04 ±0.33 ±0.29 was not measured 
in the electronic mode. 

Table 6 compares the values of the form factor ratios ry and r2 measured by E791 
in the electron, muon and combined modes with previous experimental results. 
The size of the data sample and the decay channel are listed for each case. All 
experimental results are consistent within errors. The comparison between the 
E791 combined values of the form factor ratios ry and r2 and other experimental 
results is also shown in Fig. 16 together with theoretical predictions. 

TABLE 6. Comparison of E791 results with previous results. 



Exp. Events ry = F(0) /Ai (0) ra = A2 (0)Mi (0) 

E791 6000 (e + At) 1.87 ± 0.08 ± 0.07 0.73 ± 0.06 ± 0.08 

E791 3000 (/z) 1.84 ±0.11 ±0.09 0.75 ± 0.08 ± 0.09 

E791 3000 (e) 1.90 ± 0.11 ± 0.09 0.71 ± 0.08 ± 0.09 

E687 [35] 900 (At) 1.74 ± 0.27 ± 0.28 0.78 ± 0.18 ± 0.10 

E653 [36] 300 (At) 2.00l[j |^ ± 0.16 0.82t°-22 ± 0.11 

E691 [37] 200 (e) 2.0 ± 0.6 ± 0.3 0.0 ± 0.5 ± 0.2 



The FOCUS experiments anticipate measuring r2 and ry to better precision than 
previous experiments [34] 

CHARM LIFETIMES 

The study of the charm particle lifetimes is motivated by two main goals: to 
extract partial decay rates and to study decay dynamics. The total decay width 
can be expressed as a sum of the three possible classes of decays, so the lifetime of 
a particle can be written as 

r = A^ ^_ (7) 

Ttot ^leptonic ± ^ SL ± ^nonleptonic 

The leptonic partial width is normally very small {Tieptonic ~ 10^^ — 10^^) due to 
helicity suppression. The observed TsiiD^) and Tsl{D^) are equal to within 10%, 



as expected from isospin invariance. So, the large difference observed in the D'^ 
and hfetimes, t{D^)/t{D^) = 2.55 ± 0.04, is due to a large difference in the 
hadronic decay rates (Tnonieptonic) for the and the D^. If there were no other 
diagram but the spectator and no QCD effects we would have the same value for 
r. Thus, the different lifetimes are an indication that we need to take into account 
contributions from diagrams where the W interact with two valence quarks, such as 
W-annihilation (WA) and W-exchange (WX), and any interference between them, 
as shows the Fig. 17. 
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FIGURE 17. Hadronic decays diagrams for charm meson decays. 

A systematic approach now exists for the treatment of inclusive decays, based 
on QCD and consists of an Operator Product Expansion (OPE) in the Heavy 
Quark mass. In this approach the interaction is factorized into three parts: weak 
interaction between quarks, perturbative QCD corrections and non-perturbative 
QCD effects. The decay rate is given by 



HQ 



^^S/,|V^Q,J [A, + — + — + ...] 



where A is taken as the heavy quark mass and fi is a phase space factor. Ai = 1 
contains the spectator diagram contribution; A2 is the spin interaction of the heavy 
quark with light quark degrees of freedom inside the hadron and produces differ- 
ences between the baryon and meson lifetimes; ^3 includes the non-spectator W- 
annihilation, IV-exchange and Pauli Interference (PI) of the decay and the spectator 
quarks contributions. 



New results on charm lifetime measurements are shown in table 7. The most 
relevant information is the D'^ lifetime from E791 [38,39] and FOCUS [40] which 
is now conclusively above the lifetime, with a ratio 



Rr = t{D+)/t{D^) = 1.22 ± 0.02 



(9) 



It is important to note that Rr is now ten standard deviations away from unity, 
indicating that although not dominant, the WA/WX contributions are significant. 
In fact the OPE model predicts Rr = 1.00 - 1.07 without WA/WX contributions, 
allowing a variation of ±20% if the WA/WX operators are included, in agreement 
with the new result. 

In the baryon sector SELEX and FOCUS have preliminary measurements of 
the A^ lifetime, as shown in table. 7 and Fig. 18. Using 100% of their data in 
A+ pK-n+, SELEX find r = 177 ± 10{stat)fs [41], in a 2 - 3a disagreement 
with PDG98 and FOCUS. 




FIGURE 18. D+ (left) and A+ (right) lifetimes from FOCUS (preliminary results) 



TABLE 7. Summary of new charm lifetime measurements. 

Experiment r(£)+) fs r(£>") fs T{Df ) fs t(A+) fs 

PDG'98 1057 ±15 415 ± 4 467 ± 17 206 ± 12 

E791'' 1065 ±48 413±3±4 518 ±14 ±7 

CLEO 1033.6 ± 22. 1+^2% 408±4.lti^-^ 486.3 ± 15.0l:5 i 

FOCUS'' 506±8 ' 204.5 ±3.4 

SELEX'^ 177 ± 10 

World Average 1052 ± 12 412.8 ± 2.7 499.9 ±6.1 201.9 ±3.1 

'^r(_D+) using only the 07r+ mode. 

^ Preliminary results with no systematic uncertainty quoted 



Lifetime differences and D^D^ mixing 

E791 has published searches for a hfetime difference between the CP — even and 
CP — odd eigenstates of the [39]. To do so they compared the hfetimes of the 
decays K-K+ {CP = +1) and D° ^ K-7i+, (CP mixed) [42], shown in 

Fig. 19. 




FIGURE 19. Invariant mass for D° K-tt+ (a) and D" R- K+ (b) decays (left). Ex- 
ponential fits for number of decays as a function of reduced proper decay time (right), (E791 
results). 



Defining 

T{K-K+) - TiK-n^ 



^(i^-7r^ 



ycp (10) 



The time integrated ratio of mixed to non mixed decay rates in charm meson is 
given by 

P _ nD'^D~'^f) _ x^ + y' .... 

where 



_ Am _ Ar 



with 

Am ^ mi -1712, Ar = ri-r2, r = (ri + r2)/2 

where Fi is for CP even states and r2 for CP-odd states. 
Fi apphes to D° — > K~K^ and V apphes to K~7r'^ if CP is conserved. 

Mixing can appear if there is either a difference in the masses of the CP eigenstates 
Am or if there is a difference in the decay rates AF. 

E791 observed no difference in hfetimes and quoted T{K7r) = 0.413±0.003±0.004 ps 
t{KK) = 0.410±0.011±0.006 ps, ycp = 0.008±0.029±0.010 or -0.04 < ycp < 0.06 
( 90% CL) or AF= 2{rKK - ^Kn) ^ (0.04 ± 0.14 ± 0.05) ps'^ 
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